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ABSTRACT 

We present the results of a ground-based search for the secondary eclipse of the 
3.3 Mj up transiting planet CoRoT-2b. We performed near infrared photometry using the LIRIS 
instrument on the 4.2 m William Herschel Telescope, in the H and K s filters. We monitored the 
star around two expected secondary eclipses in two nights under very good observing conditions. 
For the depth of the secondary eclipse we find in H-hauA a 3 a upper limit of 0.17%, whereas we 
detected a tentative eclipse with a depth of 0.16±0.09% in the X s -band. These depths can be 
translated into brightness temperatures of Tjj <2250 K and = 1890^5° K, which indicate 
an inefficient re-distribution of the incident stellar flux from the pla net's dayside to its nightside. 
Our results are in agreement with the CoRoT opt ical measurement ( Alonso et al. 20091 ) and with 
Spitzer 4.5 and 8 jxm results ([Gillon et al.ll2009ch . 



Subject headings: planetary systems - techniques: photometric 



1. Introduction 

The observation of a secondary eclipse, or oc- 
cultation of an exoplanet by its host star, is 
one of the few techniques that allows to ex- 
tract information about the surface or atmo- 
spheric conditions of these bodies. In the in- 
frared, it has already been possible to detect the 
thermal emission of several Hot Jupiters (start 
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2005), and reveal some features of their atmo- 
spheres, such as thermal inversions at high at- 
mospheric altitudes in planets th at receive very 



strong flux es from their host star ([Knutson et al 



2008L l2009h . These studies were only possible us 



ing the Spitzer and HST space telescopes. Using 
new observing techniques or new instrumenta- 
tions, the precisions of ground-based photomet- 
ric time-serie s have recently shown a significant 
improvement ( Alonso et al.l 2008ri Johnson et al 



120091 : iGillon et alJl2009aHSouthworth et al.ll2009l) . 

Very recently, the Na absorption on one exoplanet 



was detected from the ground (Rcdficld et al 



2008), and three independent teams announced 
ground-based detections of the secondary eclipses 
of d ifferent exoplanets (ISing &: Lopez-Morales 



20091: Ide Mooii fc Snellenll2009tlGillon et alJ2009bll . 

in z and K s bands, thus extending the studies of 
exoplanets' atmospheres to ground facilities. The 
near infrared windows close to 2 /im allow to ex- 
plore a wavelength domain that falls between the 
detections in the optical a ccessible now to sp ace 
mission s such as CoRoT ( Baglin et al.l 120021 ) or 



Kepler ([Borucki et al.l 2003) , and the wavelengths 
above 3.5 iim in th e reach of Spitzer, for which 
Gillon et al.1 (|2009cf ) recently reported measure- 
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ments on CoRoT-2b. Thus, high precision near 
infrared measurements have the potential to pro- 
vide additional constraints to the models of the 
atmospheres of these objects. 

Orbiting an active G star with a 1.74 d 
peri od, the massive pla net (3.3Mj up ) CoRoT- 
2b (jAlonso et al.l l2008al ) shows an intriguingly 



large radius of 1.47i?j up . While several of the 
Hot Jupiters show this anomalous and currently 
poorly understood bloating, CoRoT-2b poses ad- 
ditional challenges to most of the proposed mech- 
anisms to provide extra sources of energy (e.g 



Bodenheimer et al.l200ll:IShowman fc Guillotteooij : 



Winn fc Holmanl 120051 : Burrows et al.l 120031 ) due 
to its large mass. For instance, invoking increased 
opacities or kinetic energy transp ort in the plane- 
tary interior fails for CoRoT-2b (jGuillot fc Havell 
20091 ). According to these authors, the models 
point towards a young age for the planet /star 
system, and tidal dissipation due to either circu- 
larization or synchronization of the planet. The 
effect of on-going tidal circularization might be 
revealed as a shift of the secondary eclipse time 
with respect to its expected position according to 
a circular orbit. 

Due to its short period and its large radius, 
CoRoT-2b is among the exoplanets with the 
largest expected secondary eclipse depth, and 
should thus be detectable from ground based fa- 
cilities. This fact, and the unsolved issues with 
the radius of the planet, motivated us to perform 
the observations described in this work. 

After describing the way the observations were 
performed in the next Section, we present the data 
analysis in detail in Section 3, which led to an 
upper limit for the secondary eclipse depth in H- 
filter and a tentative detection of the secondary 
eclipse in K s . Section 4 will be devoted to the 
discussion of our results and their implications. 

2. Observations 

We observed CoRoT-2b during the two nights 
of 12th and 19th June 2008, using the LIRIS near- 
IR im ager and spectrograph ( Acosta-Pulido et"al~l 
20021 ) on the 4.2 m William Herschel Telescope on 
La Palma island. The observations were scheduled 
around the new moon, and the weather conditions 
were very good, with a seeing varying smoothly 
between 0.45" and 0.95" on the night of June 




Fig. 1. — A sample raw image of the observed 
field in the K s filter, showing the location of the 
target, the reference star used, and the two check 
stars that were used in Sect. 13.21 to estimate the 
effect of the systematic noises. 
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Fig. 2. — Light curves of CoRoT-2b and the refer- 
ence star in the H filter, before correcting for the 
changing airmass (top, black points) and after cor- 
rection (bottom, grey points). Vertical lines show 
the expected moments of eclipse ingress, center 
and egress. See Section |3~T1 

12th, and between 0.6" and 1.6" on June 19th. 
We used the image mode of the LIRIS instrument, 
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with a K s filter ( 1.99-2. 3 1/im) on the first night, 
and a H filter (1.47— 1.78/zm) on the second one. 
In order to avoid saturation of the detector and 
to minimize the effect of detector cosmetics, we 
used the same technique as described for obser- 
va tions with th e 1.5 m Carlos Sanchez Telescope 
in lAlonso et al. ( 2008bl ). that is, we severely defo- 
cussed the telescope, and obtained a ring-shaped 
PSF with an outer radius of about 24 pixels (6"). 
The exposure time was 15 s on the first night and 
17 s on the second one, with a typical overhead to 
read the detector and store the data of about 3 s. 
The flux on individual pixels of the target's PSF 
was never above 15K ADUs, a dynamical regime 
where the detector is known to have a good lin- 
earity. Therefore, we avoided the necessity to per- 
for m non-linearity corrections of the data as done 
by Ide Mooii k Snellen! (|2009h . In total, we ob- 
tained 1380 useful frames on the night of June 
12th, and 1096 on June 19th, observing in both 
nights for about 7.5 h around the expected times 
of secondary eclipses. A sample K s filter image is 
presented in Figure [TJ showing the position of the 
target and the reference and check stars that were 
used during the analysis. 

Following again lAlonso et aT l (|2008bh . we chose 
not to use any dithering pattern, but instead took 
twice each night a series of 10 images with the star 
shifted by 30" (offset images), in order to moni- 
tor the evolution of the image background, keep- 
ing the same exposure time. We also obtained de- 
tector calibration frames following standard tech- 
niques in IR photometry: before sunset, series of 
30 images of the dome with the dome-light 'on' 
and 'off' were used to build a master flat-field im- 
age, by subtracting a combination of the images 
taken with the lights 'on' to the combination of 
the images with the lights 'off'. Unfortunately, we 
obtained noisier results when we attempted to cor- 
rect our images by a division of this flat field, and 
we thus chose not to perform flat-fielding correc- 
tions. 

3. Analysis 

In order to search for the signal of the secondary 
eclipse in the t wo data sets, we adap ted the IDL 



routines used in lAlonso et al.l (|2008b| ). The offset 
images of each night were median-combined. We 
subtracted the mean sky background level from 
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Fig. 3. — Fully corrected curve of CoRoT-2b in the 
H filter. Vertical lines show the expected moments 
of eclipse ingress, center and egress. The solid line 
shows the best-fit trapezoid to the second-half of 
the secondary eclipse (points in black) . The data 
points before the center of the eclipse (grey dia- 
monds) are suspicious of uncorrected instrumen- 
tal effects and have not been considered for the 
eclipse depth determination. See full description 
in Section |3~T1 
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Fig. 4. — The \ 2 space for different eclipse center 
times and depths in the H-fdter dataset showed 
in Figure [31 indicating the 1,2 and 3cr confidence 
limits, using the third method described in Sec- 
tion [34] We conservatively interpret these results 
as a 3-cr upper limit of 0.17% in this wavelength. 
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this combined offset image, after measuring it in a 
region without stars. The remaining image /offset 
will thus account partially for the hot and dead 
pixels in the image, as well as for the effects from 
uncorrected cosmetics in the image. We obtained 
the best results (lowest dispersion in the final light 
curve) when we subtracted / ff se t from each sci- 
ence image, without performing any other stan- 
dard calibration. Cosmic ray impacts and remain- 
ing hot pixels were detected by applying the IDL's 
implementation of the top-hat morphological op- 
erator (e.g. Serra 1983), which works by apply- 
ing an opening (dilation) operation to the image, 
and then subtracting the result from the original 
image. This way, small features with high flux 
levels in the image can be identified; they were re- 
placed by a linear interpolation to adjacent pixels. 
Finally, from each science image we subtracted a 
sky background level that had been measured in 
a region without apparent stars. 
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Fig. 5. — Similar to Figure [H giving the light 
curves for the target and reference stars in the K s 
filter, before (top, black dots) and after correct- 
ing for the changing airmass (bottom, grey dots). 
Two check stars at similar distances to the refer- 
ence were also measured to validate the analysis 
method, and to measure the stability of the refer- 
ence star. See Section [3~2l 



The centroid of the target star was calculated 
by fitting a Gaussian with a non-zero revolution 
axis to the ring-shaped PSF. Since the target has 
a close companion at 4" distance with overlapping 



PSFs (1.72 mags fainter in K 8 and 1.8 mags fainter 
in H, from the 2MASS catalogue) that causes the 
centroid not to be well centered on the target, we 
added systematic shifts (of -0.7 and 2.7 pixels in 
the x and y coordinates of the detector) to the 
calculated target centroid. We also obtained pho- 
tometry from the ring-shaped psf 's of several ref- 
erence stars with similar fluxes as the target. To 
avoid the introduction of additional noise from in- 
dividual fits to the centers of the reference stars, 
we applied fixed position offsets between the tar- 
get and reference star apertures. The flux was 
then summed inside a circular aperture centered 
on each star. We tested apertures between 24 
and 32 pixels in radius, and obtained the best re- 
sults (less dispersion in the final curve) when using 
32 pixels for the H filter data and 26 pixels for the 
K s filter night, due to the slightly different defocus 
applied on the two observing nights. The residual 
background levels were measured in a ring around 
each star (between 2 and 2.5 times the aperture 
radius). The mean level of the target's flux was 
7.8 xlO 6 e~ on the night of the K s filter observa- 
tions, and 1.5xl0 7 e~ on the night of the H filter 
data. We recorded for each data point the residual 
sky background levels measured in the ring around 
each star, the initially subtracted sky background 
level, the x and y positions of the stars on the 
detector, and the airmass, in order to be able to 
study the possible correlations of the flux curves 
against those parameters. In the following anal- 
ysis though, only the values of the positions of 
the stars on the detector were used to justify the 
exclusion of the first data points obtained in the 
//-filter data. 

3.1. Analysis of the //-filter light curve 

We corrected the effect of the airmass by fit- 
ting a cubic polynomial to the off-eclipse sections 
of the light curve until the star passed the merid- 
ian around 3h UT (or 0.62 HJD in Figure [2]), and 
applied the fitted polynomial to the whole curve. 
The same procedure was used for the reference 
star, which was chosen among 8 different stars as 
the one whose curve had a similar shape as the tar- 
get. We noticed that even in the relatively small 
field of view of the instrument, the stars that were 
further from the target exhibited different behav- 
ior as the target, making them useless for build- 
ing a reference star from a combination of several 



4 



stars. The initial and corrected fluxes of the target 
and the reference star are plotted in Figure [2] 
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Fig. 6.— Light curve of CoRoT-2b in the K s filter 
around the phase of secondary eclipse. The verti- 
cal lines show the expected times of first contact, 
eclipse cente r and fourth contact, according to the 
ephemeris in lAlonso et al.l (|2008af) . The solid line 
plots the best fit model when the eclipse center 
time is fixed to ephemeris, and the dashed line the 
best fit model when the center time is also left free 
(superposed in this plot). See Section [3~2l 



Next, we corrected for the flux of the close com- 
panion to the target star, which provided an extra 
19% to the target flux in the ii-band. After the 
division of the target's flux by the reference star, 
there still remained some low frequency noise at 
a level of few 0.001 that we corrected by fitting a 
4 tft- -order polynomial to the parts outside of the 
predicted secondary eclipse. The final light curve 
is presented in Figure [3] Looking at Figure [2J the 
dispersion of the data is clearly higher from the 
beginning of the observations until approximately 
the mid-eclipse moment, coinciding with the ob- 
servations at bigger airmass (starting from 1.8) 
and more jitter in the guiding of the telescope, 
which is most notable in the y target-positions in 
our data (the standard deviation of the recorded 
y- values improves from 0.9 pixels at the beginning 
of the night to 0.3 pixels at the end). We sus- 
pect that this bigger jitter is due to tensions in 
the telescope mount when pointing at high air- 
mass. We thus discarded these points in the fur- 
ther analysis, because they are suspicious of be- 
ing affected by uncorrected instrumental or atmo- 
spheric effects (grey points in Figure [3]). Their 



inclusion in the analysis would result in a deeper 
and temporally more severely shifted eclipse than 
the one presented here. 

We applied three different methods to evaluate 
the depth and significance of the secondary eclipse 
signal. In the first method, we: 

• removed the best fit model from the dat^l, 
where the only fitted parameter was the 
depth, whereas the total duration and 
ingress/egress duration were fixed to the 
values of the primary eclipse, 

• shifted circularly the residuals: in order to 
keep the red noise of the data unaffected, 
all the time stamps ti were re-assigned to 
data points located at a randomly selected 
distance R in number of points ii+_R; when 
i + R > N, then the re-assigned time stamp 
was ti+R~N, where N is the total number of 
points, 

• reinserted the best fit model, 

• re-evaluated the fitted depth. 

After several hundred realizations (with different 
randomly selected values of R), the final depth, 
which takes into account the effect of red noise in 
the data, is 0.06±0.03%. We plot this model as 
the solid line in Figure 

The second method consisted of the fitting of 
two Gaussians to 1) the flux-distribution of points 
inside the full eclipse (138 points in the ii-filtcr 
data between HJD - 2454637 = 0.553 and 0.585, 
thus avoiding the moments of expected egress time 
in the case of a circular orbit) and 2) a subset of 
the points outside the eclipse (there are a total of 
532 points) with the same number of points as in 
1), and comparing their centers. This fit was per- 
formed to 500 subsets of 2), with randomly chosen 
starting points and widths (ranging from 0.0001 to 
0.0007 in normalized flux values) of the histogram 
bins used in the construction of the flux distribu- 
tion. The resulting depth is of 0.02±0.06%. 

Finally, we explored the x 2 distribution of resid- 
uals from a fit to a model in a grid of eclipse cen- 
ters (from -20 to +20 minutes from the expected 



2 A simplified model with only depth, center, total duration 
and ingress/egress duration of the occultation as parame- 
ters. 
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center-time) and depths (from 0.01% to 0.2%) of 
the secondary eclipse. The minimum \ 2 an d the 
1, 2 and 3-cr confidence levels are presented in Fig- 
ure HI and the best fitted depth using this method 
is 0.07±0.05%. None of the three methods pro- 
vided a statistically significant detection of the 
secondary eclipse in the H -filter, and we can only 
safely provide a 3-cr upper limit for the secondary 
eclipse depth of 0.17%. 

3.2. Analysis of the K s filter light curve 
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Fig. 7. — The \ 2 space for different centers and 
depths of the secondary eclipse in the K s -filter 
(Figure [6]), and the 1,2,3,4,5c confidence limits. 
A secondary eclipse with a depth of 0.16±0.05% 
is tentatively detected at the expected time for the 
case of a circular orbit. See Section l3~2l 
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Fig. 8. — The same analysis of Figure [7] performed 
on the two check stars. No signal is detected with 
3-cr upper limits of 0.06 and 0.05%, thus reassuring 
the detection presented in the previous Figure. 

We analyzed the K s -filter light curve in a very 
similar way as the ii-filter one. We chose an aper- 
ture with a radius of 26 pix, and we again cor- 
rected for the flux of the close companion to the 
target star, which provided an extra 16% of the 
flux. The plots before and after correction for the 



airmass with a low order polynomial (in this case 
a line provided good enough results) are shown 
in Figure [5j A careful inspection of that Figure 
reveals that there are slight differences in the be- 
haviors of the stars, for instance just after 0.67, 
that depend on the apparent distance in the sky 
among the different stars. In order to verify that 
the signal we are searching for is not due to this 
residual systematic noise, we repeated the same 
analysis with two check stars located at different 
parts of the detector and at opposite directions 
from the reference star. To reduce the high fre- 
quency noise of the reference star, we applied a 
4-points median filter. The final light curve (after 
division by the flux of the filtered reference star) 
shows a tentative secondary eclipse event (Fig- 
ure [6]). We evaluated the depth and significance of 
the secondary eclipse using the same three meth- 
ods described in the ii-filter section. Model fitting 
and shifting of residuals (method 1) resulted in a 
measured depth of 0.17±0.09%. The fit to differ- 
ent gaussians (method 2, with 246 points consid- 
ered to be completely inside eclipse: 0.56 <HJD - 
2454630< 0.61, and 801 points in the out of tran- 
sit sample) gave a depth of 0.15±0.15%, while the 
X 2 mapping (plotted in Figure [7]) resulted in an 
eclipse with a depth of 0.16±0.05%. The differ- 
ent error bars reveal the presence of systematic 
noise in the data that are intrinsically taken into 
account in methods 1) and 2), but not in the \ 2 
mapping. We chose as a final value 0.16±0.09% 
for the eclipse depth measurement, i.e., the mean 
value of the three methods. The \ 2 mapping was 
also performed on the two check stars, and the Fig- 
ure[8]shows 3-cr upper limits of 0.06 and 0.05% re- 
spectively for a secondary eclipse, re-assuring the 
detection of the secondary eclipse in the .?r s -band. 

4. Discussion 

Our analysis of the near infrared observations 
obtained at the H and K s bands has allowed to 
establish a 3-cr upper limit to the occultation of 
CoRoT-2b in the ii-band of 0.17%, while the K s 
filter observations reveal a secondary eclipse of 
0.16±0.09%. 

In order to convert the detected depths (or up- 
per limits) in the two band passes into brightness 
temperatures, we to ok a modeled spectrum of a 
G8V star from the IPicklesI (|l998l ) library. We 
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calibrated the model in order to have the same 
total flux as a 5625 K Planck function. Un- 
der the assumption of a black-body planetary 
emission, and zero Bond albedo, we computed 
the flux of the planet for different temperatures. 
Using the H and K s filter response function^, 
we could then estimate the secondary eclipse 
depth for different temperatures. The 3-ct up- 
per limit on the iJ-filter corresponds to an upper 
limit of the brightness temperature Tjj <2250 K, 
while the Jf s -band measurement implies a black- 
body emission of Tj<-,=1890l^ K. If we as- 
sume that the planet is in thermal equilibrium, 
then the measured re-distribution factor is of 
/=0.57lg'32, which points to an inefficient re- 
distribution of the incident stellar flux from the 
planet's dayside to its nightside. These results 
are in agreement, despite their bigger uncertain- 
ties, with the values obtained from the analysis 
of t he optical light curv e from the CoRoT satel- 



lite (lAlonso et all 12003) . of_Tc oRoT 



Recently, iGillon et al.l 
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1910±?° K. 



(|2009d) published 



ondary eclipse depths from Spitzer observations, 
of 0.510±0.042% and 0.41±0.11% at 4.5 and 8^m 
respectively, which led these investigators to also 
favor models of inefficient heat-distribution. 

CoRoT-2b falls in the pM-cla ss of planets 



accor ding to the classification by iFortnev et al 



(2008), which exhibit hot stratospheres (temper- 
ature inversions), caused by currently debated or 
unknown upper atmosphere absorbers. Accord- 
ing to these authors, such an effect might lead 
to the detection of "anomalously" deep secondary 
eclipses in the near infrared. In the previous in- 
terpretation, we have neglected the reflected light 
component, which is expected to be smaller than 
the thermal emission at these wavelengths. 

Additional observations will help to provide a 
firm detection in the H band, as well as to improve 
and confirm the observed K s band depth of the 
secondary eclipse. 
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